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ABSTRACT 

In this paper we consider methods of analysis and optimal design of redshift surveys. 
In the first part, we develop a formalism for analysing galaxy redshift surveys which 
are essentially two-dimensional, such as thin declination slices. The formalism is a 
power spectrum method, using spherical coordinates, allowing the distorting effects 
of galaxy peculiar velocities to be calculated to linear order on the assumption of 
statistical isotropy but without further approximation. In this paper, we calculate 
the measured two-dimensional power for a constant declination strip, widely used in 
redshift surveys. We present a likelihood method for estimating the three-dimensional 
real-space power spectrum and the redshift distortion simultaneously, and show that 
for thin surveys of reasonable depth, the large-scale 3D power cannot be measured 
with high accuracy. The redshift distortion may be estimated successfully, and with 
higher accuracy if the 3D power spectrum can be measured independently, for example 
from a large-scale sky-projected catalogue. 

In the second part, we show how a 3D survey design can be optimized to measure 
the power spectrum, considering whether areal coverage is more important than depth, 
and whether the survey should be sampled sparsely or not. We show quite generally 
that width is better than depth, and show how the optimal sparse-sampling fraction 
/ depends on the power P to be measured. For a Schechter luminosity function, a 
simple optimization fP ~ 500ft. -3 Mpc 3 is found. 



1 INTRODUCTION 



The measurement of fluctuation power in the galaxy dis- 
tribution is an important test of galaxy formation models, 
since the fluctuation spectrum, of mass at least, is predicted 
readily by such models. Power can be measured from three- 
dimensional redshift surveys, or from projected catalogues, 
by numerical inversion techniques. Redshift surveys have the 
advantage (and disadvantage) that they are distorted by the 
effects of peculiar velocities, and can therefore be used to 
extract information on the density parameter, under the as- 
sumption that structure grows by gravitational instability. 
One has then the possibility of measuring three-dimensional 
power and the density parameter (via /3 = Q,q' 6 /b, where b 
is the bias parameter for the survey in hand) from a galaxy 
redshift survey (Heavens & Taylor 1995; hereafter HT95; see 
also Kaiser 1987, Hamilton 1992, Cole et al. 1994). 

It is clear that the longest wavelength which can be 
measured is limited by the size of the survey, so it is attrac- 
tive to consider surveys which are essentially one- or two- 
dimensional, to maximise at least one dimension without 
incurring prohibitive cost in observation time (e.g. Broad- 
hurst et al 1990). The difficulty with such an approach as 
a method for measuring the power spectrum is that a low- 
dimensional power measurement at a given wavenumber will 
have a contribution (which may be dominant) from much 
smaller scales in three dimensions (e.g. Kaiser & Peacock 
1991). The interpretation of the observed power spectrum 
can therefore be difficult. For surveys which do not corre- 
spond to the 'distant-observer' approximation (cf. Kaiser 



1987), the power spectrum measurement and the redshift 
distortion become linked, and this further complicates the 
analysis. 

The ease with which the parameters of interest may be 
extracted depends on the choice of coordinate system and 
basis functions in which the density field is expanded. For 
two reasons the choice of spherical polar coordinates is com- 
pelling. Firstly the survey is almost certain to be defined in 
terms of a fixed areal coverage (independent of depth), and 
secondly, a flux limited survey will have a selection function 
<f> (or, equivalently, a mean observed density po(r)) which 
is dependent on distance, but not on direction. The mean 
density of the survey is then separable in spherical coordi- 
nates p(r) = po(r)M(0, cp), where M is either 1 or depend- 
ing on whether the direction (8, if) is in the survey or not. 
The second reason is that, unless /3 -C 1, it is impossible 
to ignore redshift distortion effects, and since the distortion 
between the real-space map and the redshift-space map is 
purely radial, it is straightforward to include the distortion 
in a power-spectrum analysis in spherical coordinates (cf 
Zaroubi & Hoffman 1996). 

The choice of basis functions must also be done with 
some care. It is very useful to choose functions which pick 
up a narrow range of wavenumbers from three-dimensional 
space. In this regard, spherical Bessel functions are ideal, 
and they also have advantages in that the redshift distor- 
tion is relatively simple in this system. An arbitrary choice 
of functions (or equivalently, an arbitrary choice of projec- 
tion onto the sky) makes the interpretation of measured 2D 
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power difficult, since the 3D wavenumbers contributing to 
the power may be poorly constrained. Even in this idealized 
case, the range of 3D power contributing to the 2D coeffi- 
cients can be large, especially for thin surveys of a few de- 
grees thickness. In these cases, 3D power spectrum at large 
scales becomes difficult to achieve with high accuracy. 

One of the aims of the paper up to this point is to 
demonstrate that the expected accuracy of a proposed sur- 
vey in measuring some parameter can be estimated in ad- 
vance, and this sort of study can and should influence how a 
survey is designed, whether in 2D or 3D. We investigate in 
section 4 how 3D design may be optimised for power spec- 
trum estimation, subject to constraints. We find that, if ob- 
serving time is constrained, then it is always better to cover 
a large area on the sky, rather than going deep. In some cases 
it can be advantageous to sample galaxies sparsely, with the 
sparse-sample fraction being dependent on the power to be 
measured and the luminosity function of the galaxies. We 
present a very simple formula for calculating the optimal 
sparse-sample fraction. 

The paper is laid out as follows: in Section 2 we con- 
sider power measurements in idealised cases of ID and 2D 
surveys, which highlight the problems which such surveys 
have to address. In Section 3 we present a new method for 
analysing constant declination strips, for measuring the real- 
space power spectrum and redshift distortion. In Section 5, 
we solve the problem of designing surveys optimised for 3D 
power estimation. 



2 ID AND 2D: GENERAL CONSIDERATIONS 

The fractional overdensity is <5(r) = p(r)/p — 1, where p is 
the mean density. Our Fourier transform convention is = 
J d 3 rS(r)e- A r with inverse 8{r) = — / d 3 k <S k e lk r . 
The Power spectrum is defined by 

(<5k<5k'} = (27r) 3 P3£i(fc) <5 D (k-k') 

so that the correlation function is 



(1) 



€(r) 



d 3 kP3D(k) e*' r . 



(2) 



(2tt) 3 

In an idealized pencil-beam survey, the density field along 
a line is measured, and the ID power spectrum estimated. 
To relate this to the 3D power spectrum, we follow Lums- 
den, Heavens & Peacock (1987), noting that the correlation 
function is the same along the line as in 3D, by isotropy. 



Pi D (k) 



j dx£(x)exp( — ikx) 
1 



(2^)3 



dxd^Pzoik) exp(ik x x) exp(-ikx) 



(3) 

where we assume the pencil beam lies along y = z = 0. The 
integration over x gives a ID delta function, and changing 
the remaining integration over k x and k y to a polar integra- 
tion, we get 

= / dkP 3D (k)k (4) 

( 27r ) J w 

Hence we see that power in one dimension comes from all 
shorter wavelengths in 3D. This is readily understandable by 
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Figure 1. Kernel function G(k, k)/k for a thin sheet (solid) and 
pencil-beam survey (dashed). The 2D wavenumber is unity. 



consideration of the following 2D — * ID illustration: imagine 
looking along a corrugated roof at an angle to the corruga- 
tions. The separation of peaks along the line is longer in 
ID by a geometrical factor, so the ID power (averaged over 
angles) has contributions from 2D power at all shorter wave- 
lengths. 

Notice also that the ID power spectrum must be a 
monotonically decreasing function of k (the squared ampli- 
tudes of the Fourier coefficients may not be monotonic, being 
drawn from a Rayleigh distribution). Also note that if the 
3D power spectrum has a cutoff at some large wavelength, 
the ID power spectrum will be constant (and non-zero) on 
all larger scales. 

For comparison with later analysis in this paper, we 
define the kernel G(k, k) such that the measured power is 



P(k) 



f 

Jo 



d\nkP 3D {k) G(k,k) 



(5) 



from which we see that the kernel for a ID skewer is 
G(k,k) =2Trk 2 G(k-k) 



(6) 

where O is the Heaviside function. This unpleasant convo- 
lution function is shown in Figure 1 (throughout we plot 
G(k,k)/k, since we use a linear rather than logarithmic k 
axis). For a practical pencil-beam survey, the kernel will be 
suppressed at high k, but this calculation illustrates the se- 
vere problems in interpreting the power spectrum of pencil- 
beam surveys - the ID power may be coming from much 
larger wavenumbers in 3D. 

For thin, infinite plane surveys, a similar analysis (Pea- 
cock 1991) gives 



P 2D (k) 



Jk 



dkP 3D (k) 



k 



^k 2 - k 2 ' 



(7) 
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The associated kernel G(k, k) = 2k 2 0(k - k)/\/k 2 - k 2 is 
also shown in Figure 1. Again what one finds is that the 
kernel feeds a lot of power (if it exists) from high k to low k. 
If the plane survey is of finite thickness, the high-fc part of 
the kernel will be suppressed, and weighting of the data can 
help still further (see Section 3) but it is still possible that 
2D power at k may have little connection with 3D power on 
such scales. 



3 2D SURVEYS 

For surveys which have one dimension considerably smaller 
than the other two, it is sensible to reduce the dimension- 
ality of the survey by projection onto a 2D surface. We can 
then reduce the dimensionality of the transform correspond- 
ingly. As emphasized in the introduction, there are consider- 
able advantages in using spherical coordinates for the anal- 
ysis. The survey is almost certainly characterized by a ra- 
dial selection function, and an angular selection defining the 
boundary (cf Fabbri & Natale 1990, Scharf et al. 1992, Scharf 
& Lahav 1993, Fisher et al 1994a,b 1995). Also, the effects 
of redshift distortion enter only in the radial direction. In 
the case of a thin slice survey, one avoids any difficulties 
which might be apparent in 'flattening-out' the survey into 
a plane suitable for cartesian analysis. One final point is that 
it is the largest scale modes where the radial nature of the 
distortion is most important to treat correctly. 

In this section, we develop a 2D expansion of the density 
field projected onto a fixed declination, allowing for redshift 
distortions. We use coordinates s,0,<p, where s = cz/Hq is 
the distance assigned on the basis of the redshift z, assuming 
uniform expansion (Ho is the Hubble constant). This is the 
normal distance assigned in redshift surveys, and it differs 
from the true distance r because of peculiar velocities. In 
this paper we consider only z«l, but the spherical coordi- 
nate system allow the effects of non-euclidean geometry and 
temporal evolution to be included if desired. 

Our 2D expansion is based on the 3D Fourier-Bessel 
expansion of the density field (cf HT95, Lahav 1993): 



Stm{k) 



with inverse 



d i v8{v)j l {kr)Y l m *{9, V ) 



5(r) = 



\J\ Z)/ dkk2s (m(k)j e (kr)Y e m (9,ip). 



(8) 



(9) 



The statistical properties of 8i m (k) are derived in the ap- 
pendix. We proceed by projecting all galaxies onto the cen- 
tral value 9 — 9o, and expand in terms of m and k. At this 
stage, we leave the choice of radial expansion function gen- 
eral, f{kr). We also allow for radial and angular weighting 
of the data via the functions w s (s) and wn(Q), which may 
help in optimizing the signal-to-noise and apodizing. The 
radial weight may be fc-dependent. For a constant declina- 
tion strip, the obvious orientation of coordinates is to have 
the centre of the strip at constant 9 — 9o, and to expand in 
terms of m. We let the thickness be A9, and the width Aip, 
centred on tp = 0. Our choice of expansion is 



Pm{k) 



— I d sp(s) f(ks) exp(— irrup) w s (s)wci(£i). (10) 



As an important aside, there is an issue over which frame of 
reference should be used for redshift-space expansions of this 
sort. Should the redshift be measured in the Local Group 
frame or the Microwave Background frame? In either case, 
the redshift distance is 



s(r) = r 



(v - v ) • r 



H r 2 



(11) 



where vo is the peculiar velocity of the frame of reference. 
This relationship is general, but since we wish to make a per- 
turbation expansion, we must ensure that the second term 
in the square brackets is always small. Assuming a suffi- 
ciently coherent velocity field such that v approaches the 
Local Group velocity as r — > 0, we see that the expansion 
must be done in the Local Group frame. 

The difference between the expansion coefficients and 
their mean values 



Pm(fc) 



2 / 

d rpo(r) f(kr) exp(— imtp) w 3 (r)wn(fl) (12) 

7T ' 



can be related to the Se m {k) by substituting for S(r) from (B), 
and noting that number conservation implies that p(r)d 3 r = 
p(s)d 3 s: 



D m (k) 
where 



p m (k) - p m (k) 



dk Si m (k) Aj(fc, k) k 2 



(13) 



(21+1) (I-m)\ f ,(*+|*|) /a 
4tt (£ + m)P ' 



2 sin [(m - m)A</p/2] 
(m — m) 



cos(e -Ae/2) 



cos(0 o +Ae/2) 



tfcifV) (14) 



and 

A|(fc, k) 
*/(*,*) 

Vi{k,k) 



*f(fc,fe)+/3V*(fc,fc) 
2 

TV 



dr po(r) ji(kr)f(kr) w s (r) r 



2 
nk 2 
d_ 
dr 



d 

drp (r)— [f(kr)w s (r)] x 
dr 



j- [j- e (kr)] r 2 . (15) 
The signal part of the covariance matrix can be written as 

(D m (k)D* m ,(k')) = J dlnkP 3D (k)G mm ,(k,k',k) (16) 

where 

G mm , (k, k',k) = (2vr) 3 zp m 'Ai(k, k)Ae(k\ k)P (17) 



T - WP™WP m *. The shot noise contribution 
to the covariance matrix is 



and Zf m 



(D m {k)D* m ,{k')) S N = - J drd9d<t> 

r 2 sin 8po(r)f(kr)f(k'r) exp \i{m — m')(p\ w 2 (r)wa(Q,) . 

(18) 

In practice, one splits D m (k) into real and imaginary parts, 
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with similar, but more cumbersome, expressions for the co- 
variance matrix elements. Some kernels are shown in Fig. 
2 and 3 for a survey with a gaussian selection function 
exp [— (r/r*) 2 ] , with r* = 450/iT 1 Mpc, with survey limits 
A6> = 6° centred at declination 30°, and width Acp = 90°. 
The radial expansion function is chosen to be a spherical 
Bessel function, with 1 = 2. This choice is motivated in two 
ways. Firstly, we know that in 3D the Bessel functions give 
narrow kernels, so they seem a good start in 2D. A sec- 
ond, related point, is that the function gives little weight to 
the very nearby part of the survey. Since this is the thinnest 
part, it is likely to contribute significantly to aliasing difficul- 
ties. The weighting scheme chosen also helps in this regard. 
Fig. 3 shows kernels for plane waves with almost the same 
wavenumbers as Fig. 2, and with direction along the central 
ip value, to correspond as closely as possible to Fig. 2 (cf the 
analysis of the Las Campanas survey by Landy et al. 1996). 
These curves are simply integrals of the squared modulus of 
the window function transform, calculated via a 200 3 FFT, 
which accounts for their slightly ragged nature. The com- 
parison between the methods is not quite straightforward, 
as the 2D modes look rather different. Note that the Fourier 
modes assume f3 = 0; for non-zero /3, the kernels are ex- 
tremely complicated in the Fourier case. The comparison is 
most stark if one compares the dimensionality of the objects 
which one needs to calculate to include redshift distortions. 
In essence, the 2D coefficients are linear combinations of the 
3D coefficients. If no simplification is possible, one needs 
to calculate a 5D object to calculate a range of 2D coeffi- 
cients. This is required if one uses Fourier modes (Zaroubi & 
Hoffman 1996, equation 10), but using the spherical modes 
reduces the dimensionality such that the most complicated 
objects are only 3D (see (13)). It is this fact that the kernels 
for non-zero j3 can be readily calculated for the spherical 
modes which is their major advantage. It arises, of course, 
from the radial nature of the distortion and selection func- 
tion, and the use of angular coordinates to delimit the sur- 
vey. 

For high m the spherical kernels are sometimes not cen- 
tred on k, and the 2D power may come principally from 
shorter wavelengths in 3D. The effects of this, and shot 
noise and cosmic variance can be accounted for correctly 
using likelihood techniques, so the 3D power spectrum can 
be estimated, but it is clear from Figs. 2 and 3 that the 
task is not going to be easy, whichever method is used. The 
accuracy with which the power and j3 determination can be 
done with the Fourier-Bessel transform is explored in the 
next sections. 

3.1 Parameter estimation 

We can use the analysis method presented in the last section 
to estimate the real space power spectrum and (3 maximizing 
the likelihood. Symbolically 

W,P) = (2n) N ^\\C\\ 1 / 2 6XP \-\Y, D ^A ■ ( 19 ) 

where C is the covariance matrix of the N data values, de- 
pendent on the (parametrized) P(k) and /3. Once again, the 
data are the real and imaginary parts of D m (k). The like- 
lihood method has the advantage that all the aliasing ef- 
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Figure 2. Kernel function for constant declination slice, 
for modes with m = 2, and, from left to right, k = 
0.008, 0.016, 0.025/1 Mpc" 1 . Solid lines arc for (3 = 1, dotted 
for /3 = 0. 
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Figure 3. Kernel function for Fourier modes with same wavenum- 
bers as Fig. 2. 
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fects are treated correctly, and we do recover 3D power esti- 
mates with correct error bars. We illustrate this method by 
analysing a numerical simulation created with Couchman's 
AP 3 M code (Couchman 1991). A slice between declinations 
20 and 40 degrees is projected onto declination 30, with a 
right ascension range of 90 degrees. The power spectrum is 
a power-law P(k) oc k~ , and the likelihood for the am- 
plitude of P(k) and (3 is shown in Figure 4. The details 
of the analysis are that the nonlinear wavenumber (where 
k 3 P(k) I (2tv 2 ) = 1) is 183 (units are arbitrary) , and the anal- 
ysis examines modes up to fcmax = 30. Pushing the maxi- 
mum analysed wavenumber beyond this pushes (3 down, as 
the effects of Fingers-of-God become apparent. These could 
be reduced by including the effects of smoothing (cf. HT95) , 
but they have not been incorporated here. The m modes are 
analysed in steps of 2 from 2 to 20, and the wavenumbers 
selected are from 6 to 30 in steps of 2. There is no difficulty 
in principle in taking every m and n mode, but there is a 
numerical problem as adjacent modes are too strongly cor- 
related and the matrix becomes numerically singular. Note 
that modes separated by 2 are correlated, and the corre- 
lations correctly accounted for in (19). We also put a con- 
straint on the wavenumber perpendicular to the line-of sight, 
ensuring that it is not too nonlinear, by rejecting modes with 
mfc/4.0 < fcmax- This ensures that the transverse wavenum- 
ber at the peak of j2{kr) is no more than 4.0/3.3 times fc max . 
Experimentation shows that this gives unbiased estimation 
of (3 and P(k). The galaxies are weighted with the Feldman 
et al. (1994) optimized weighting w s (r) = [1 + po(r)P(k)]~ 
and P(k) in the weighting is taken as a constant, compara- 
ble to the true power in the simulation. The true parameters 
are shown by the encircled cross. We see that the method 
is capable of determining f3 and the power spectrum with 
somewhat larger errors than a fully 3D survey (HT95) with 
similar numbers of objects, and also note that here we need 
to examine larger wavelengths than in the 3D case (up to the 
nonlinear wavenumber/6, as opposed to nonlinear wavenum- 
ber/3 in the 3D case). Failure to do this leads to underesti- 
mation of f3 because of nonlinear effects. 

3.2 Errors on f3 and P(k) 

In this section, we calculate the expected errors in a deep, 
thin-slice survey, such as might be achievable in the first year 
of the AAT 2dF survey. The error on the parameters is read- 
ily estimated using the Fisher information matrix (Tegmark, 
Taylor & Heavens 1996). 

For a set of parameters 6i, i — 1,N (e.g. (3 and the 
power spectrum in N — 1 wavenumber bins), the covariance 
matrix of the parameter estimates is 

T = (66 T ) - (6)(6 T ) = F -1 (20) 
where 

Fa = ^Tracc (C _1 C,iC _1 C,j) (21) 

is the Fisher information matrix. C is the covariance matrix 
of the data (DD T ), and C ; = dC/dOi. The Fisher matrix 
is readily obtained from the data covariance matrix (16). To 
illustrate this, we calculate the parameter covariance matrix 
for a thin slice, 6° x 90° , with a Gaussian selection function 
po(r) = p„ exp(-r 2 /r»). We take p* = 0.02h 3 Mpc -3 and 
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Figure 4. Likelihood function for the amplitude of the real-space 
power spectrum and j3 for a numerical simulation whose true pa- 
rameters are shown by the cross. The contours are separated by 
0.5 in ln(likclihood). 



r» = 450ft -1 Mpc, broadly comparable to an optical survey 
to a limit b — 19.5 (cf forthcoming AAT and Sloan galaxy 
surveys) . 

We analyse modes from m = 2 to m = 20, once again 
separated by 2 to avoid the covariance matrix becoming 
numerically singular. The k values are spaced by 0.0167 h 
Mpc -1 , and the modes are analyzed up to k = 0.05/i Mpc -1 , 
consistent with our previous numerical experiments for un- 
biased results. The summations extend to I = 60, and the k 
integrations extend to k = 0.165/i Mpc -1 . The galaxies are 
weighted with the Feldman et al. (1994) optimized weight- 
ing w s (r) = [1 + po(r)P(k)]~ and P(k) in the weighting is 
taken to be 2700 h~ 3 Mpc 3 . The expected error on (3 from 
such a slice is 0.236, and the expected fractional error in 
the power is shown in Fig. 5, for a power spectrum assumed 
to be smooth on a scale of 0.0167 h Mpc -1 . Increasing the 
width of these k bins decreases the error. Note how the error 
increases at the high-fc end beyond the maximum wavenum- 
ber analysed (4.0/3.3 times 0.05 h Mpc -1 ~ 0.06), and at 
the low-end, where the size of the survey becomes compa- 
rable to the wavelength (2iv/r r ~ 0.014). The correlation 
matrix for the parameters is shown in Table 1. This anal- 
ysis takes only a matter of minutes on a workstation, once 
the matrices $ and V have been calculated. These take a 
few hours, but are calculated once only for a given survey. 
What is apparent from this example is that, even for a deep 
survey with many objects, P(k) is detected on scales of the 
survey k ~ 27r/450 ~ 0.014fo Mpc -1 , but not with good ac- 
curacy. j3 estimation is actually not bad (error 24%), but 
this could be improved noticeably (to 15%) if the 3D power 
spectrum is determined independently, from a sky-projected 
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Figure 5. Expected fractional error on P(k) from a deep thin 
slice survey, for narrow bins in fe-space. For details, see text. 

catalogue such as the APM survey. An application of this 
method to the Las Campanas survey (Shectman et al. 1995) 
is in progress. 
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Table 1. Correlation matrix for the parameters, in the or- 
der P and the five fractional power spectrum measurements 
in order of increasing k. 



4 3D SURVEYS: OPTIMIZING FOR POWER 
ESTIMATION 

We have shown in the previous section how the expected 
error on a parameter may be estimated in advance for a 
given analysis method and survey design. The conclusion 
that 2D surveys are not particularly good for determining 
large-scale 3D power suggests that genuine 3D surveys may 
be more profitable. The issues of design and analysis are also 
relevant in 3D, and here we consider the problem of opti- 
mising the design of a galaxy redshift survey to measure the 
power spectrum on some particular scale. The typical deci- 
sions to be made are whether to go for a deep survey over 



a small area of sky, or a shallower survey over a wider area 
of sky. We also consider whether it makes sense to sparse- 
sample the galaxies, or to observe every one. This section is 
essentially a Fourier analogue of Kaiser's (1986) treatment of 
sparse-sampling to estimate the two-point correlation func- 
tion, generalized to account for a radial selection function, 
and with the proper power error estimate of Feldman, Kaiser 
& Peacock (1994) incorporated. 

For simplicity, we make the following assumptions: we 
assume that the effects of redshift-space distortions arc 
small, and we assume that the observing time for each galaxy 
is proportional to the inverse square of its flux. The former 
is motivated by earlier studies (HT95) where the optimal 
weighting was found to be insensitive to the degree of red- 
shift distortion (see also Hamilton 1996). The latter assump- 
tion is an example; different constraints, for fibre systems for 
example, could be incorporated if desired. Our prime con- 
straint is that the total duration of observing is taken to be 
fixed. In considering sparse-sampling, we restrict attention 
to a sparse-sample fraction which is constant for all galaxies 
in the parent sample. Thus, for example, we do not consider 
a variable sparse-sampling rate which depends on flux. 

It was shown by Feldman et al. (1994), HT95 and 
Hamilton (1996) that the optimal weighting of galaxies in 
the survey is 



w (r) 



l + /n(r)P(fc) 



(22) 



where P(k) is the (prior estimate of) the power to be mea- 
sured, and n(r) is the mean number density of galaxies at 
position r. We introduce the possibility of sparse-sampling 
by multiplying this number density by a factor /. 

Feldman et al. demonstrated that this weighting gives 
rise to an error in the power of cr 2 P /P 2 = (2-7r) 3 /(Vk-0, where 
Vk is the volume of k— space over which the power is aver- 
aged, and 



I(f,S,Q)=Q / dr 



( 1 + /pi(r)) 



(23) 



Here S and Q are the flux limit and solid angle of the survey. 
Our problem then reduces to maximising I with respect to 
/, S and fl, subject to the constraint that the total observing 
time is fixed. 

To do this optimisation, we need the luminosity function 
$(!/), from which the number density is obtained: 



n(r, S)=n (X) 



dL $(L). 



(24) 



where X = 4irr 2 S. If the time to observe an object of flux 
density S' is A/(16tt 2 ,S" 2 ) for some constant A, then the time 
to observe a fraction / of all objects to a flux limit S in a 
solid angle f2 reduces to 



t(/,5,n) = Qf / drXr 
Jo 
QXf 



poo />c 

/ dr \r 

JO J 4tt 



dL 



dXX 5/2 n 2 (X) 



£(L) 
L 2 



2(4nS)?/ 2 J 



where n 2 (X) = dL $(L) L~ 



(25) 
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The time constraint then simply yields S oc (fi/) 2//r , 
and the error is minimised when 



p/7 



dX- 



x 1 ' 2 



i 1 + /Pno(X)) 



(26) 



is maximised. f2 and / may be chosen freely, apart from 
the obvious limits, with the depth of the survey S being 
dependent on the choice. We see immediately that the er- 
ror is minimised if O is made as large as possible. This is a 
quite general result, consistent with the general knowledge 
that surveys should be wide before being deep. If we fix the 
solid angle of the survey (as large as convenient), then we 
can straightforwardly solve for / to optimize the error. The 
analysis is readily generalized for observing times which are 
proportional to S'~ a r~ r3 (a = f3 = 2 might be appropriate 
for fixed- width slit spectroscopy). In this case, one maxi- 
mizes 



' l + 2a/3-/3/3 



/ 



dX- 



x l/2 



1 



fPn Q (X) 



(27) 



Fig. 6 illustrates the effect for a Schechter luminos- 
ity function $(L)dL = 0* (L/L*)" 1 ' 3 exp(-L/L*)dL/i*, 



with <f>* = 0.013/1 3 Mpc -3 . 



The optimal sampling occurs at 



fP ~ 500/i~ 3 Mpc 3 , although of course / itself is bounded 
above by unity. To the left of the minimum, shot noise be- 
comes dominant, whereas to the right, the extra sampling 
reduces the volume observable, so that cosmic variance dom- 
inates. To estimate /, the power at a wavenumber k = 0.01 
to O.lh Mpc -1 is about 1000-10000 h~ 3 Mpc 3 , depending on 
the galaxy type and theoretical prejudice (e.g. Baugh & Efs- 
tathiou 1993,1994, Ballinger et al. 1995), which motivates a 
sparse-sampling strategy of / ~ 0.1. The error rises rapidly 
if fP < 100, so one must take care not to under-sample. If 
the power spectrum on large scales has the Zel'dovich form 
P oc k, a survey to measure very large-scale power should 
be sampled fully. 



5 CONCLUSIONS 

In this paper we have presented a new method for analysing 
thin, near-constant declination slice surveys, using a 2D pro- 
jection and expansion in radial and angular functions. We 
have also considered the optimisation problem of depth and 
sparse-sampling for 3D surveys. There are two main advan- 
tages of using spherical coordinates for analysis. The first is 
that the survey is usually defined by a fixed areal coverage, 
and a flux limit which leads to a selection function which is 
purely radial. The second advantage is that the effect of red- 
shift distortion is radial, so it is straightforward to include 
it in the analysis. By expanding in spherical and angular 
functions, one can treat linear redshift distortions without 
further approximations, and this allows, in particular, anal- 
ysis of long-wavelength modes which do not subtend small 
angles on the sky. By using carefully chosen radial functions 
for the analysis, one can ensure that the modes one analyses 
essentially include only 3D modes which are still linear. 

For an all-sky survey, the formalism leads to very sim- 
ple analysis, and this method is clearly the best we have to 
date. What was not clear was whether the method could 
be adapted for surveys of relatively small areal coverage, 



Power spectrum variance 




fP/(h 



3„„ i(; 3\ 



Figure 6. The unnormaliscd variance in the power spectrum for 
a fixed solid angle and a Schechter luminosity function, as a func- 
tion of the sparse-sample fraction / and the power spectrum to 
be estimated. The details of the luminosity function are given in 
the text, and the units of P are /i -3 Mpc 3 . 



since the mixing of modes of different I and m make the 
expansion more cumbersome. This paper shows that, even 
with thin, essentially 2D surveys, one can retain the advan- 
tages of the spherical expansion without severe additional 
complexity. Our error analysis shows that 3D power can be 
estimated from 2D surveys, properly including the effects 
of aliasing, mode-mode correlations, shot noise and cosmic 
variance. However, the reduction in dimensionality means 
the errors achievable are unlikely to be very small. This is in 
contrast to 3D surveys, where small errors on the real-space 
power spectrum can be achieved from the Fourier-Bessel 
technique (Ballinger, Heavens & Taylor 1995). In 2D, /3 can 
be determined with reasonable accuracy (about 25%), but 
the best approach will probably be to use a sky-projected 
catalogue to estimate the 3D power independently, and then 
to use spherical harmonics with the slice to measure (3 with 
higher accuracy (about 15%). 

We also show in this paper that 3D surveys may be op- 
timised for measuring 3D power, given a constraint on total 
observing time, by choosing as wide an area of sky as pos- 
sible, and by sparse-sampling at a rate which is dependent 
on the expected power to be measured. 



APPENDIX 

In this appendix, we calculate the covariance matrix for the 
continuous spherical transform. We transform the density 
field <5(r) = ^— ^3 J d 3 k<5 k e lkr and expand the exponential 
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as a sum of Bessel functions 
8{r) = j dkdV.uh > 



J2i e Mkr)Y e m (e k ,^)Yr(9,^k 2 



(28) 



where Q k = (6k,fk) is the direction of the k vector, and 
(6, if) is the direction of r. The definition of the spherical 
harmonics used in this paper is that found in Binney & 
Tremaine (1987): 



47r (£+|m|)! 



x cxp(im0) x 



(-1)' 



m > 
m < 



(29) 



The spherical expansion is 
^2 



km{k) = \j- J S(r)j e (kr)Y e m (9, l p)d i r (30) 
which becomes 

8 lm {k) = (2n)-' i/2 i e JdkdQ l ,S k Yr(e k ,^)S D (k-k) (31) 

where we have used the orthogonality (Binney & Quinn 
1991) 



/ 



dQdrje (k'r)je(kr) r 2 Y e 7'*(6, <p) Y e m (9, <p) 
-^g8 D (k - k ') s te> s mm> 



(32) 



where 8 D and 8 K are Dirac and Kronecker delta functions 
respectively. The covariance matrix is 

(8 lm {k)8l, m ,{k')) = (2n)- [i i 1 - 1 ' J dkdk'dSl- k dn- k , 

{5&,)Yr (e- k ,^)Yir'(9- k „v- k ')8 D (k - k)8 D (k' - k'). 



(33) 



Defining the power spectrum by 
<«> = (2tv) 3 P(k)S D (k-k') 

S D (k-k') !;D 



(2*yp{k) 



fc 2 



- ^w)8 D {vv -fi'), 
(34) 

(n = cos 6) the analagous expression of orthogonality for the 
Fourier-Bessel modes is 

(5 em (k)5p m ,(k')) = P{k) SD{k k ~ k,) 5^,8^,. (35) 

The power is evenly divided between real and imaginary 
parts, except for m = modes, which are real. 

Redshift distortions 



In a redshift space map, galaxies are placed at a position 
s—(s,6,(f>), where the distance coordinate s is the recession 
velocity divided by the Hubble constant Ho- In general this is 
not the true distance because the galaxy may have a peculiar 
velocity v. The redshift space position is then related to the 
real-space position r by 



s(r) 



1 + 



(v - vo) ■ r 



(36) 



where vo is the peculiar velocity locally. 

To expand the spherical expansion to linear order, we 
first note that p(r)d 3 r = p(s)d 3 s, and make a Taylor expan- 
sion of the resulting integrand to first order in s — r 

j e (ks)w s (s) ~ j e (kr)w s (r) + (s - r)-^- \je(kr)w s (r)] . (37) 

dr 

To obtain an expression for s — r, we assume potential flow 
v = — V$ (valid for linear, growing- mode perturbations), 
where 4>(r) is the velocity potential. The effect of the local 
group velocity vo is to add an extra term to the mean of 
the transform coefficients, and will be treated separately. 
Expanding $ in terms of $e m (k), we find 



s — r = 



v.r 

So 

1 



2 E 

Cm 



dk$> lm {k)^±Yr*{e^)k 2 . 



(38) 

The peculiar Poisson equation V 2 $ = /35(r) relates the 
potential to the galaxy overdensity field, which leads to 
&im(k) = —f38 lm (k)/k 2 . From this we find (choosing units 
such that Ho = 1) 



(3 



TV 



dk8 em (k)^^Yr{e, V ). (39) 
dr 



This leads to the V matrix terms in the main text. The 
effect of the local group velocity is to add the following to 
the mean value of D m (k): 



w s (r)f(kr) 

dQ. wn(£i) exp(— imip)\o ■ f 
where r and vo are unit vectors. 



(40) 
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